Incommensurate magnetism in cuprate materials 
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In the low doping region an incommensurate magnetic phase is observed in La2-xSr x CuOi. By 
means of the composite operator method we show that the single-band 2D Hubbard model describes 
the experimental situation. In the higher doping region, where experiments are not available, the 
incommensurability is depressed owing to the van Hove singularity near the Fermi level. A pro- 
portionality between the incommensurability amplitude and the critical temperature is predicted, 
suggesting a close relation between superconductivity and incommensurate magnetism. 



PACS numbers: 74.72.-h, 75.10.-b, 71.10.Fd 

The dynamical spin susceptibility for cuprate mate- 
rials has been investigated by inelastic neutron scatter- 
ing and NMR techniques. Neutron scattering data on 
La,2-x{Ba, Sr) x CuOi have shown [1-7] that away from 
half-filling the commensurate antiferromagnetic phase is 
suppressed and short-range incommensurate antiferro- 
magnetism develops. The magnetic Bragg peak in the 
dynamical structure factor S(k, ui) broadens and devel- 
ops a structure with four peaks located at [(1 ± S)tt,tt] 
and [tt, (HS)tt]. The incommensurability amplitude 8{x) 
does not depend on the frequency and temperature; it is 
zero up to the doping x w 0.05 where a commensurate- 
incommensurate transition is observed, then, increases 
with the hole concentration x, following the linear law 
8(x) ~ 2x up to x Ri 0.12; beyond this point there is 
a deviation downwards [7]. Unfortunately, experimental 
data above x = 0.18 are not available, due to the dif- 
ficulty in preparing single crystals. It is important to 
stress [5] that the value of doping x = 0.05, where the 
transition is observed, corresponds to the value of doping 
where the material becomes superconducting. These in- 
commensurate spin fluctuations are not observed in other 
cuprate materials; a flat topped magnetic peak has been 
observed [8] in Y Ba-iCu-iO^+y with y « 0.6, while for 
the case of the electron-doped Nd,2- x Ce x Cu04 no in- 
commensurate magnetism has been observed [9]. The 
difference in the spatial modulation experimentally ob- 
served in La2~ x Sr x CuOi and Y Ba^Cu^O^+y has been 
related to a difference in the topology of the Fermi surface 
[10,11]. 

From a theoretical side the presence of incommensu- 
rate magnetic correlations was firstly found in the 2D 
Hubbard [12] and in the t — J model [13] by Quantum 
Monte Carlo (QMC) calculations. To improve the sit- 
uation the t — t' Hubbard model has been considered, 
but the results are not definite and there is no general 
agreement [14]. 

In this Letter we shall study the spin magnetic sus- 



ceptibility of the two-dimensional single-band Hubbard 
model by means of the Composite Operator Method 
(COM) [15-16]. We have shown [16] that in the static 
approximation the dynamical spin susceptibility is given 
by the following expression: 
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where n is the particle density, D is the double occu- 
pancy and by Q^^k, ui) we mean the retarded part of 

Qa/J 7 *(k»w) =i J ^G Q/3 (k+p,^ + r!)G 75 (p,r!). The 
2x2 matrix G(k, ui) is the thermal causal Green's func- 
tion, defined by G(k, ui) — {T[tp(i)ip^(j)])F.T., where ip(i) 
is the doublet composite operator 
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with = c a (i) [1 - n_o- (i)] and r/ a (i) = c CT (i)n_ CT (i). 

c a (i) is the electron operator at site i. By means of the 
equation of motion and by considering the static approx- 
imation, where finite lifetime effects are neglected, the 
Green's function G(k, u)) can be computed in the course 
of a fully self-consistent calculation where no adjustable 
parameters are considered [15]. 

The calculation of the uniform static susceptibility 
Xo(x,T) for various values of doping and temperature 
has been given in Ref. 16, where we showed that Eq. 
(1) qualitatively reproduces the experimental situation 
observed in La,2~ x Sr x CuO\. The principal results ob- 
tained in Ref. 16 can be so summarized. For a fixed 
temperature, \o(x,T) is an increasing function of the 
doping, reaches a maximum at a critical doping x c , then 
decreases. The value of x c does not change with tempera- 
ture and is determined by the ratio U /t, varying from to 
1/3 when U/t changes from zero to infinite; for U/t = 4 



1 



we have x c — 0.27. This doping dependence qualita- 
tively reproduces the experimental behavior observed in 
La 2 - x Sr x CuOi [17], where a critical value x c « 0.25 is 
reported. For a given doping, when the system is away 
from the critical density, xo{ x -,T) as a function of tem- 
perature has a behavior similar to a 2D Heisenberg an- 
tifcrromagnet with a maximum at a certain temperature 
T m ; the position of T m decreases as the system is doped 
away from half-filling and tends to zero for x — ► x c ; in 
the vicinity of the critical doping there is a large increase 
of Xo{x, T) for low temperatures. The behaviors of T m as 
a function of x and of Xo(x, T) as a function of T m well 
reproduce the experimental data of Ref. 17. 

The peak exhibited by Xo(x, T) for a certain critical 
doping is related to the fact that for x — x c the Fermi 
energy crosses the vHs. This can be seen in Fig. 1, where 
N(Ep), the density of states calculated at the Fermi en- 
ergy, and Xo(x,T), the uniform static susceptibility, are 
given versus the doping parameter x. We have chosen 
U/t = 4 and k B T/t = 0.01. At half-filling the Fermi 
energy is at the center of the two Hubbard bands; by 
varying the dopant concentration some weight is trans- 
ferred from the upper to the lower band, Ep moves to 
lower energies and crosses the vHs for a critical value of 
the doping; further increasing x, Ep moves away from 
the vHs. A study of the Fermi surface shows that for 
x > x c we have a closed surface which becomes nested 
at x = x c and opens for x < x c . An enlarged Fermi 
surface with a volume larger than the noninteracting one 
has been reported by QMC calculations [18,19] and by 
other theoretical works [20] . 

To understand the role played by the vHs in the case of 
spin fluctuations, it is useful at first to consider the case 
of noninteracting Hubbard model (i.e. U = 0). What 
we learn [21] from the study of this model can be so 
summarized. 

- The k-dependent susceptibility x(k) exhibits a maxi- 
mum at a certain value k* which depends on the position 
of the vH energy u v h with respect to the Fermi energy 
E F . 

- When n — 1 u) v h = Ep and k* = Q = (n, tt). In 
addition there is a singularity coming from the nesting of 
the Fermi surface and the staggered susceptibility x(Q) 
exhibits a stronger divergence than xo- 

- When n / 1 the Bragg peak at Q = (tt, tt) opens in 
four peaks, situated at k* = [tt(1 ± <5) , tt(1 ± 5)]; there 
is a transition from commensurate to incommensurate 
magnetism. The incommensurability amplitude S(x) in- 
creases as a function of the doping x with the same law 
as the shifting of Ep with respect to the vHs: 



uj v h - Ep w ax 4 / 3 
S(x) w bx^ 3 



a « 3.62 
b w 1.31 
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the susceptibility. Calculations based on the use of Eq. 
(1) show that around Q = (tt, tt) x(k) has an incommen- 
surate structure along the four corners of a square, with 
a minimum at Q . This incommensurate structure con- 
tains a mixing of two components. The relative position 
and the intensity of the two contributions change signifi- 
cantly with doping. To study the effect of the interaction 
we see that the /c-dependent susceptibility x(k) can be 
written as x(k) = Y%,j=i Xij( k ) where 



X«(k) 
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Kij(Kp) (4) 



In the interacting case, the shifting of the vHs and the 
band structure have a drastic influence on the form of 



The quantities iWj(k, p) are expressed in terms of the 
spectral intensities. The term Xinter — X12 + X21 de- 
scribes transitions between the two bands £i(k), the up- 
per Hubbard band, and ^(k), the lower Hubbard band; 
while the two terms X11 and X22 describe intraband tran- 
sitions. Since E\ takes values mostly above the chemical 
potential, the contribution of X11 is small. The interband 
term Xinter is reported in Fig. 2. This term originates a 
peak in the susceptibility, which moves from the commen- 
surate position Q = (tt, tt) to (tt, tt/2) when the doping is 
increased from x — to x = x c . The intensity of the peak 
decreases by increasing doping. The intraband term X22 
is reported in Fig. 3. This term gives a peak which is a 
reminiscent of the Van Hove singularity in the density of 
states. At zero doping the vHs is far from the Fermi en- 
ergy and the peak is located (tt/tt/2) and has a low inten- 
sity. When doping increases, the peak increases its inten- 
sity and moves along the line (k x — tt,tt/2 < k y < 3tt/2). 
At the critical doping x — x c the vHs lies on the Fermi 
energy and the Fermi surface is nested. Then, the peak 
°f X22 is situated at Q and has a very high intensity, due 
to the concomitance of these two effects. It is interesting 
to note that the peak position of X22 moves towards Q 
with the same law as given in Eq. 3. 

The total susceptibility is reported in Fig. 4 for three 
values of doping. For zero doping we mainly have a com- 
mensurate structure with a peak coming from Xinter , lo- 
cated at (tt, tt), and a smaller peak, coming from X22, lo- 
cated near (tt,tt/2). Upon doping, the two peaks moves 
for different reasons. X22 moves because the Van Hove 
singularity moves towards the Fermi energy. Xinter moves 
because the band structure changes with doping. When 
the critical doping X — X £ IS reached, the Van Hove sin- 
gularity is at the Fermi energy and the Fermi surface is 
nested. A commensurate structure is recovered with a 
very high peak coming from X22 • 

In Fig. 5 the incommensurability amplitude 5(x) is re- 
ported as a function of doping. In the region of low (high) 
doping the peak coming from Xinter (X22 ) is predominant 
and very well separated from the other; in these regions 
5(x) has been evaluated as the middle point of the half- 
width of the peak. In the region 0.10 < x < 0.18 the two 
peaks overlap and 5(x) has been calculated by taking the 
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average of both peaks and we have a plateau due to the 
superimposition of X12 and Xinter- For comparison we 
report the experimental data of Refs. 4, 5 and 7. The 
linear behavior of S(x), observed in the low doping region, 
agrees exceptionally well with the experimental data, re- 
ported in Refs. 4-7; the downward deviation reported in 
Rcf. 7 for x > 0.12 might correspond to the plateau the- 
oretically observed. One of the most striking feature of 
the results presented in Fig. 5 is the similarity between 
the incommensurability amplitude S(x) and the critical 
temperature T c . 5(x) is maximum in the region of opti- 
mal doping where T c is maximum. It has already been 
observed in Ref. 7 that there is a linear relation between 
S(x) and T c up to the optimal doping level x ~ 0.15. 
Our theoretical results extend to the all region of doping 
a relation of proportionality between S(x) and T c . 

The same result for 5(x) can be obtained by consid- 
ering Im x(k, u)). Some results have been given in Ref. 
22. We preferred to study the k-dependent susceptibility 
x(k) because this quantity provides more strict informa- 
tion about the spatial range of the magnetic correlations. 
On the other hand an exact experimental determination 
of x(k) is not easy, since it must be calculated by the 
accessible S*(k, u) through a Kramers-Kronig relation. 

The present analysis shows that the interaction in the 
Hubbard model has mainly two effects. One is the change 
of the critical doping from x = to some critical x c , 
due to the shift of the vHs. This shift explains and well 
reproduces the unusual normal state behavior of xo m 
hole-doped cuprates. The other is a band structure effect 
which is responsible of the incommensurate modulation 
of the magnetic susceptibility in the low doping region. 

The picture that emerges is that the magnetism probed 
by neutron scattering experiments is correlated with the 
carrier density. In the low doping region the suscepti- 
bility is mainly controlled by the term Xinter which de- 
scribes band structure effects and then reflects the topol- 
ogy of the Fermi surface. In the overdoped region the 
Fermi energy is close to the vHs and the effect of nesting 
in the intraband term is important. In Y Ba^Cu^O^+y 
we have a different topology of the Fermi surface and no 
nesting is expected; this might be the reason why incom- 
mensurability is not observed. 

The main results obtained in this Letter can be so 
summarized. There is experimental evidence that in 
hole-doped high T c cuprates the Fermi level is close 
to the vHs for values of doping close to those where 
the superconducting phase is suppressed. In the con- 
text of the Hubbard model a van Hove scenario well 
describes some of the unusual properties observed in 
the normal state, but an analysis show that this sce- 
nario is related to the overdoped region and not to 
the optimal doping. The existence of a critical doping 
where the vHs lies on the Fermi energy should imply 
a peak in the staggered susceptibility. Then, we pre- 
dict that commensurate magnetism should be recovered 



in the nearness of the critical doping, implying a propor- 
tionality relation between the incommensurability ampli- 
tude S(x) and the superconducting critical temperature. 
Recalling that in Lai- x Sr x CuO^ the commensurate- 
incommensurate transition is observed at the same value 
of doping x ~ 0.05 where superconductivity starts, at 
least for La,2- x Sr x CuOi, a scenario [23] which relates 
the superconducting phase to the presence of incommen- 
surate magnetism emerges. 
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Figure Captions 

Fig. 1 The density of states at the Fermi energy N(Ep) 
and the uniform static susceptibility Xo( x i t) as functions 
of the doping x. U/t = 4 and k B T/t = 0.01. 

Fig. 2 The interband term X3(k) along the line k = 
(■K,k y ) for ksT/t = 0.01 and for various values of the 
doping x < 0.27 with step 0.03. U/t = 4. 

Fig. 3 The intraband term X2(k) along the line k = 
(it, k y ) for ksT/t = 0.01 and for various values of the 
doping x < 0.27 with step 0.03. U/t = 4. 

Fig. 4 The spin magnetic susceptibility x(k) along 
the line k = (n,k y ) for various values of the doping x. 
U/t = 4 and k B T/t = 0.01 . 

Fig. 5 The incommensurability amplitude 5{x) as a 
function of the doping x. The dashed line indicates the 
theoretical result. U/t = 4 and k B T/t = 0.01. 
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